PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



1 



Reconstruction of IceCube coincident events and study of 
composition-sensitive observables using both the surface and deep 

detector 

Tom Feusels*, Jonathan Eisch^ and Chen Xu*, for the IceCube Collaboration^ 

*Dept. of Subatomic and Radiation Physics, University of Gent, B-9000 Gent, Belgium 
^Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 
^Bartol Research Institute, University of Delaware, Newark, DE 19716, USA 
§See the special section in these proceedings 



Abstract. The combined information from cosmic 
ray air showers that trigger both the surface and 
underground parts of the IceCube Neutrino Observa- 
tory allows the reconstruction of both the energy and 
mass of the primary particle through the knee region 
of the energy spectrum and above. The properties 
of high-energy muon bundles, created early in the 
formation of extensive air showers and capable of 
penetrating deep into the ice, are related to the 
primary energy and composition. 

New methods for reconstructing the direction and 
composition-sensitive properties of muon bundles are 
shown. Based on a likelihood minimization pro- 
cedure using IceCube signals, and accounting for 
photon propagation, ice properties, and the energy 
loss processes of muons in ice, the muon bundle 
energy loss is reconstructed. The results of the high- 
energy muon bundle reconstruction in the deep ice 
and the reconstruction of the lateral distribution of 
low energy particles in the surface detector can be 
combined to study primary composition and energy. 
The performance and composition sensitivity for both 
simulated and experimental data are discussed. 
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I. Introduction 

The cosmic ray spectrum covers many orders of mag- 
nitude in both energy and flux. In the energy range ac- 
cessible to the IceCube Neutrino Observatory (~0.3 PeV 
to 1 EeV) the slope of the spectrum remains mostly 
constant, except for a feature at around 3 PeV where 
the spectrum steepens. This feature is called the knee 
of the cosmic ray spectrum, and its origin is unknown. 
Proposed explanations include changes in acceleration 
mechanisms or cosmic rays leaking from the galactic 
magnetic field starting at this energy. Measurement of 
mass composition in this range could give clues to 
the origin of these cosmic rays. The IceCube detector, 
together with the IceTop air shower detector, provides 
an opportunity to measure the composition of cosmic 
ray particles in the region of the knee and beyond. 

The IceTop detector, high on the Antarctic Plain at 
an average atmospheric depth of 680g/cm 2 , consists of 



a hexagonal grid of detector stations 125 m apart. Each 
station consists of two ice tanks that act as Cherenkov 
media for measuring mainly the electromagnetic compo- 
nent of cosmic ray air showers. In each tank, two Digital 
Optical Modules (DOMs) are deployed, which contain a 
10 inch PMT and digital readout and control electronics. 
The primary energy can be reconstructed by the IceTop 
surface array (TJ. 

Deep below each IceTop station is a string of the 
IceCube detector with 60 DOMs evenly spaced between 
1.5 and 2.5 km in the ice. Combined, the IceTop and 
IceCube arrays can reconstruct the air shower core 
position and direction while measuring the shower signal 
strength at the surface and the energy deposition of 
the high-energy muon bundle in the deep ice. With 
these measurements, the energy and mass of the primary 
cosmic rays can be reconstructed. 

A characteristic difference between showers induced 
by light and heavy nuclei for a fixed primary energy 
is their number of muons, with higher mass primaries 
producing more muon-rich showers. However, the muon 
multiplicity is not directly measured by either the IceTop 
or IceCube detectors. The deep IceCube detector is 
sensitive to Cherenkov light coming mainly from energy 
loss processes of high-energy muon bundles. This energy 
loss is a convolution of the muon multiplicity of the 
shower, the muon energy distribution and the energy 
loss of a single muon. If the energy loss behavior of 
muon bundles can be reconstructed accurately, it can be 
used as a primary mass indicator (2). In Section ITVl the 
reconstructed muon bundle track described in Section Hill 
will be used as a seed to reconstruct the muon bundle 
energy loss. Simulated data is compared to experimental 
data in Section [V] to examine the detector performance 
and its sensitivity to composition. 

II. Data Sample and Simulation 

Our experimental data sample was taken from the 
month of September, 2008. At that time, the detector 
consisted of 40 IceTop stations and 40 IceCube strings. 
A total livetime of 28.47 days was obtained by selecting 
only runs where both detectors were stable. The events 
were processed at the South Pole with a filter which 
required at least three triggered IceTop stations and at 
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least 8 triggered IceCube DOMs. After this filter about 
3.31-10 6 coincident events remained. 

To study the direction and energy reconstruction of 
muon bundles in the ice, a large number of proton 
and iron showers between lOTcV and 46.5 PeV were 
simulated with the CORSIKA package. SIBYLL 
2.1 0, was used as the high-energy (>80TcV) 
hadronic interaction model, while FLUKA08 (6), Q 
was used as the low energy hadronic interaction model. 
For this study, CORSIKA was configured to use a model 
of the South Pole atmosphere typical for the month 
of July [8|. The showers were simulated according to 
an E _1 spectrum and then reweighted according to an 
E~ 2 - 7 spectrum before the knee (at 3 PeV) and an E~ 30 
spectrum after the knee. 

The IceCube software environment was used to re- 
sample each CORSIKA shower 500 times on and around 
the detector, to propagate the high-energy muons through 
the ice, and to simulate the detector response and trigger. 
The simulation was filtered the same way as data and 
yielded about 9.010 4 proton and 9.0-10 4 iron events. 

III. Direction Reconstruction 

The direction reconstruction by IceTop will be de- 
scribed first because it will be used later as a seed for 
an IceCube muon bundle reconstruction algorithm. 

An initial shower core position and direction is de- 
termined using the extracted times and charges of the 
recorded pulses from the IceTop tanks. The first guess 
core position is the calculated center of gravity of tank 
signals, while the initial direction reconstruction assumes 
a flat shower front. This shower core and direction are 
then used as a seed to fit the lateral distribution of pulses 
with the double logarithmic parabola (DLP) function 
described in [9|. Because the reconstruction of the core 
and direction are highly correlated, the resolution is 
improved by fitting the shower core position and the 
shower direction together, with a curved shower front 
and the DLP lateral particle distribution. This fit, and 
some geometrical quality cuts discussed later, gives an 
angular resolutiorQ of 1 .0° for iron showers (see the dot- 
dot-dashed line in Fig. [TJ and a core resolution of 15.0 m 
(see dot-dashed line in Fig. 0. 

The resolution depends strongly on where the shower 
core lands with respect to the IceTop array. This analysis 
uses only events with a shower core reconstructed within 
the geometrical area of the IceTop detector. The recon- 
structed muon bundle track was also required to pass 
within the instrumented volume of the IceCube detector. 

The direction that was already reconstructed by the 
IceTop algorithms alone can serve as a seed for an algo- 
rithm more specialized in reconstructing muon bundles 
in the ice. This algorithm uses only the charges measured 

'The resolution of observable Y is defined according to 
f* p(AY)dAY = 0.683, where x is the resolution and p(AY) 
the frequency of distances between the true and the reconstructed 
observable. 
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Fig. 1. The distribution of the angles between the true direction 
and the reconstructed direction for simulations of iron showers is 
shown for different algorithms. The dot-dot-dashed curve shows the 
reconstruction which uses IceTop information alone. For the dot- 
dashed line, the reconstructed core position by IceTop is fixed and the 
zenith and azimuth are determined by using a muon bundle algorithm 
seeded with the track determined by IceTop. The solid line illustrates 
the ideal limit of this reconstruction method by using the true shower 
core position instead of the reconstructed one. The dashed curve 
is obtained when a second iteration between IceTop and IceCube 
algorithms is used. 



by the IceCube DOMs and takes into account the range- 
out of the muons in a bundle 1101 . By keeping the 
reconstructed core position on the surface fixed, a large 
lever arm of at least 1500 m is obtained. This limits 
the track parameters (zenith and azimuth) during the 
minimization procedure and reduces the number of free 
track parameters from 5 to 2. 

In Fig. [1] it can be seen that this method (dot-dashed 
line) improves on the angular resolution determined by 
IceTop alone. If the core position can be determined 
more accurately, the direction reconstruction will be 
even better. The solid line on Fig. Q] represents the ideal 
limit, obtained using the true core position. Therefore, 
to improve the core resolution the new direction is kept 
fixed and used to seed the IceTop lateral distribution 
function which then only fits the core position (dashed 
line on Fig. |2). Iterating over both the surface and the 
deep detector reconstructions with progressively better 
core position and direction seeds, leads to the optimal 
resolution. The ideal limit for the core resolution is 
acquired by seeding the IceTop algorithm with the true 
direction and is shown on Fig. [2] (solid line). After the 
second iteration this limit is already obtained and gives 
a core resolution of 14.0 m and an angular resolution of 
0.9° (see dashed lines in Figures Q] and [2]). 

Using this combined method for muon bundle di- 
rection reconstruction, an almost energy independent 
angular resolution of 0.8° and a core resolution of 
12.5 m is obtained for proton induced showers. An im- 
provement of the core resolution and direction resolution 
also improves the resolution of the reconstructed shower 
size and shower age. 
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Fig. 2. The distribution of distances between the true shower core 
position and the reconstructed core position at the surface for iron 
showers is plotted for different reconstructions. For the dash-dotted 
line the IceTop DLP function that determines the core position used a 
first guess direction. The dashed line, where the IceTop algorithm was 
seeded with a better direction reconstruction from the deep ice, gives 
a slightly better core reconstruction. When the true direction is used 
as a seed, the ideal limit represented by the solid line is obtained. 



Reconstructed slant depth behavior 
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Fig. 3. An example of reconstructed muon bundle energy loss for a 
single 17 PeV iron and proton shower. The reconstructed slant depth 
behavior follows the true energy loss reasonably well. The triangles 
(squares) are the muon bundle energy loss processes for a certain 
proton (iron) 17 PeV shower with a zenith angle of 28.4° (12.6°). 
The spread of the points illustrates the stochastic nature. The solid 
(dashed) line is the reconstructed muon bundle energy loss function, 
described in the text, for the proton (iron) shower. 



IV. Muon Bundle Energy Reconstruction 

An energy reconstruction algorithm for single muons 
was previously reported in Ref. fl2j . Using IceCube 
signals and lookup tables together with a previously 
reconstructed track a constant energy loss is fit with a 
likelihood function. The lookup tables model the South 
Pole ice properties and the propagation of Cherenkov 
photons through the ice. A single high-energy muon 
above 730 GeV loses energy mainly by radiative pro- 
cesses like Bremsstrahlung and pair production, which 
produce secondary electromagnetic cascades in the ice 
along the muon bundle track. Therefore, an infinite light 
source with mono-energetic cascades every meter is used 
as a model for a single muon. 

This light model can also be used for muon bundles. 
The main difference is that a slant depth dependent 
energy loss will be needed because of the range-out of 
muons. Here, the slant depth is defined as the distance 
along the muon bundle track between the point where the 
muon enters the ice and the point where the Cherenkov 
light is emitted. 

The equation for the muon bundle energy loss is: 



fdE. 



V dX 



(X) = 



Bundle 



where 4^- is the energy distribution of the muons and 



E m ^ lt) dN dE 

(1) 



-jrjjp is the energy loss of a single muon. -Eminfsurf) = 
| [e bx — l) is the minimum energy that a muon needs 
to get to depth X. -E max ( sur f) oc ^ is the maximum 
energy a muon from a shower induced by a particle with 
A nucleons and primary energy Eq can have. 

Using a simple power-law as an approximation for the 
Elbert formula lfT3l . which describes the multiplicity of 



high-energy muons in air showers, the differential muon 
energy distribution becomes: 



k(A ) (B> 



.ST 7 "" 1 



(2) 



where 7^ = 1.757 is the muon integral spectral index 
and k is a normalization that depends on the shower 
properties. 

With the solution of the single muon energy loss 
equation, E^X) = (j5 M (surf) + § ) e" bx - f, the 
average energy loss formula can be expressed as a 
function of the muon energy at the surface: 



dX 



(X) = -a-bE^X) 



b (£7 |i (surf) + - e 



-bX 



(3) 



with a = 0.260 GeV/m, the ionization energy loss 
constant and b = 0.000357 m _1 , the stochastic energy 
loss constant from ifTTI . 

The average muon bundle energy loss function is then 
obtained by integrating Eq. ([T} using Eqs. (f2]i and (0. 
This energy loss fit function, with k and Eq/A as free 
parameters, will be used to scale the expected charges 
in a DOM from the likelihood formula in lfl2l instead 
of scaling it with a constant energy loss. 

In Fig. [3] the curves show the reconstructed muon 
bundle energy loss functions for 17 PeV primaries. The 
data points are the energy losses calculated from simula- 
tions. It can be clearly seen that the muon bundle energy 
loss function obtained by minimizing the likelihood 
function describes the depth behavior for these two 
showers better than a constant energy loss function. 

It has been shown in [2 | that proton and iron show- 
ers are separated better by the muon bundle energy 
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Comparison of Simulation and Data 
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Fig. 4. The reconstructed muon bundle energy loss evaluated at a 
slant depth of 1650 m versus the reconstructed shower primary energy. 
Shown are the center 68% of the proton simulation (right diagonal 
hashes), the iron simulation (left diagonal hashes) and the center 68% 
of the data (enclosed in the rectangles). The median of the distribution 
is also shown for the simulated data sets. The events were filtered as 
described in Section [II] with the additional quality requirement that the 
muon bundle reconstruction algorithm used signals from at least 50 
DOMs, to remove events with poorly fit energy loss. This quality cut 
removes a higher portion of muon bundles with low energy loss and 
will be accounted for in a full composition analysis. 

loss at smaller slant depths. The reconstructed IceCube 
composition-sensitive parameter which will be used fur- 
ther on in the coincidence analysis, is the energy loss 
at the top of the IceCube detector, at a slant depth of 
1650 m. At this slant depth, Cherenkov light from show- 
ers with a zenith angle up to 30° can still be detected 
by the upper DOMs, making the energy reconstruction 
more accurate over the entire zenith range. 

V. Combined Primary Mass and Energy 
Reconstruction 

Fig. |4] shows a comparison of proton and iron primary 
simulation and experimental data described in Section 
ITU using the reconstruction methods described in the 
sections [III] and [TV] While this plot has only rough 
quality cuts, it can be seen that the spread in the 
simulation is similar to the spread in the data. The 
median muon bundle energy loss from an iron primary 
shower is approximately a factor of 2 higher in the 
ice than the median muon bundle energy loss from 
a proton primary shower for the same reconstructed 
primary energy. There is a large overlap area, where 
shower to shower fluctuations overcome the effect of the 
primary mass. It is difficult to reconstruct the primary 
mass of a single shower with any certainty due to these 
large fluctuations. 

VI. Conclusion and Outlook 

The IceTop and IceCube detectors of the IceCube 
Neutrino Observatory can be used together for an im- 
proved air shower core location and direction reconstruc- 
tion. The promising method of reconstructing the muon 



bundle energy loss behavior will be further developed 
to be used in a measurement of cosmic ray mass and 
energy. This study used only a subsample of the available 
data; with more statistics and an enlarged detector these 
methods can be extended up to 1 EcV. 
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